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Abstract 

We discuss the bremsstrahlung production of soft real and virtual photons in a quark-gluon 
plasma at thermal equilibrium beyond the Hard Thermal Loop (HTL) resummation. The physics 
is controlled by the ratio Q 2 /q 2 of the virtuality to the energy. When Q 2 jq 2 <C g 2 , where g is the 
strong coupling constant, the emission rate is enhanced by a factor 1/g 2 over the HTL results due 
to light-cone singularities and the bremsstrahlung is induced by scattering of the quark via both 
transverse and longitudinal soft gluon exchanges. When Q 2 /q 2 increases, the enhancement factor 
is given by q 2 /Q 2 . When this ratio is near unity, the bremsstrahlung contribution is of the same 
order as the rate predicted by the HTL resummation. In that case, the bremsstrahlung is induced 
by both soft and hard gluon exchanges. 
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1 Introduction 



Several methods have recently been used for calculating the rate of soft real or virtual photons in a 
hot quark-gluon plasma. These methods lead to drastically different physical pictures and orders of 
magnitude for the calculated rates. In the following we will always assume the plasma to be in thermal 
equilibrium and the QCD coupling constant g <C 1 . We furthermore assume that the produced photon 
does not interact with the quark-gluon plasma due to the smallness of the electric charge e compared 
to the strong interaction coupling: in other words the photon is not thermalized. 

In the hard thermal loop (HTL) resummation analysis (l], [2|, § the photon is radiated off in processes 
which always involve an interaction between the plasma particles which are dominated by soft (energy 
~ gT) quasi- fermion exchanges P]. In the case of real photon emission, a strict application of the HTL 
rules leads to a prediction with a logarithmic singularity associated to collinear divergences [||, || . 

In a semi-classical approach |7|, ||, B], the photon is radiated off by a hard quark (energy ~ T) 
interacting with quarks and gluons in the plasma via the exchange of soft gluons, the interactions 
being assumed to be screened by the Debye mass. Multiple re-interactions of the radiating quark in 
the plasma regularize potential collinear singularities in the lowest order process and, more importantly, 
suppress the photon spectrum at small frequencies due to the Landau-Pomeranchuk-Migdal (LPM) 



effect [P-Ofl . For a different approach of this effect, see |ll|, |12|. 

The regularization of divergences in the bremsstrahlung process has also been achieved, under 
somewhat different hypotheses, by the resummation of soft photon emission with a single interaction 
of the radiating quark in the plasma |l3|, [T^| . We will not consider this possibility as it affects the 
photon spectrum at very small momenta, of order eT, while we are interested here in momenta of 
order gT. 

In the following, we come back to the problem of real or virtual soft photon emission, in the 
framework of thermal field theory, showing the necessity of going beyond the HTL approximation 
for obtaining consistent predictions. The picture which emerges is rather different from the above 
quoted results. Compared to the HTL approach we find, in fact, that there is no contribution from 
hard thermal loops at the expected leading order in g. However some non leading terms in the HTL 
approximation become enhanced by singularities near the light-cone and dominate over the would-be 
HTL order In the case of real or quasi-real photons (virtualities Q 2 /q^ <C 1) the dominant 

diagrams are the bremsstrahlung ones in agreement with the hypotheses of ||, BJ and we find that, 
for Q 2 /qQ < g 2 , the production rate is R ~ e 2 g 2 T 4 /q 2 l modulo some logarithmic factor. The results are 
shown to be independent of the covariant gauge parameter. Unlike in ||, |9| we treat the interaction 
between quarks and gluons in the plasma "exactly", in the framework of thermal field theory, and 
we do not find that the static approximation is a good one as exchanges of "magnetic" gluons are as 
important as that of "electric" ones. 

In summary, the result of our analysis concerns photons of momentum much less than T radiated 
by hard quarks in a plasma: for virtualities Q 2 /q^ <C 1 we find an enhancement, due to light- 
cone singularities, over the rate R ~ e 2 g 4 T 4: /qQ (modulo logarithmic factors) expected in the HTL 
approximation; for virtualities Q 2 jq\ ~ 1 no enhancement occurs but the bremsstrahlung contributions 
are of the same order as the "soft-fermion" contributions calculated in Q. This extends out of the 
light-cone our previous results on photon production on the light-cone (jDJ. 

Problems in the HTL expansion associated to light-cone singularities have already been discussed 
several times in the literature: as mentioned above they render the rate of soft real photon production 
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logarithmically divergent p, ^, [Tq| . In fact, already in their original paper M Braaten and Pisarski 
discuss mass divergences, i.e. terms such as g 2 T/ (to — p) with uj — > p, and they mention that these 
terms could spoil the HTL expansion. In order to solve this problem it has been proposed recently 
[jl7| , [i~8| , |l~9f to extend the HTL action to include thermal corrections (essentially the effective thermal 
mass) on hard propagators but no general proof exists that this improvement is sufficient or complete. 
A particularly interesting feature of the problem we consider here is that the mass singularities change 
the expected order of the leading HTL prediction and enhance it by a factor which may be as large 
as 1/g 2 in some cases. 

2 Discussion of the existing approaches 

This part is devoted to a critical appraisal of the previous calculations of the real and virtual photon 
production rates in a quark-gluon plasma. 

2.1 Weaknesses of the Hard Thermal Loop expansion 

One of the first applications of the HTL resummation scheme was the calculation of the production 
rate of a soft photon, of virtuality Q 2 , at rest in a quark-gluon plasma Before resummation 
is taken into account, the photon is produced by the annihilation of a soft quark- antiquark pair in 
the plasma. The imaginary part of the vacuum polarization diagram, which is proportional to the 
invariant production rate R up to a statistical weight, is found to be 

Imn%(Q) ~e 2 Q 2 (l-2n F (|)) -eV^ 2 for q,q ~gT (1) 

In the HTL approach the diagrams to be calculated are shown in Fig. |l[ In fact, only the first one, 
with effective quark propagators and effective q — q — 'j* vertices, is said to contribute: the self-energy 
diagram with the gluon loop in fact does not exist because there is no g — g — 7 effective vertex 
(Furry's theorem) and the two tadpole diagrams vanish because the traces of the corresponding 4- 
point effective vertices also vanish. A compact expression, involving convolutions of the effective soft 
fermion spectral functions has been obtained by Braaten, Pisarski and Yuan Q and one can estimate 
the order of magnitude of their result to be as in Eq. ([l]). An important fact concerns the statistical 
factor (1 — 2n F (^)) associated to the annihilating q and q which becomes (1 — n F (uj) — n F (^ — u)), 
where u is the energy of one of the annihilating soft fermion in the HTL approach: for soft energies 
it brings an extra suppression factor g to the photon production rate. 

Consider now, for example, the last diagram with the gluon tadpole and in particular the order of 
n 00 or n* J . Using the HTL rules, the 5 — 5 — 7 — 7 effective vertex is found to be of order e 2 while 
the soft effective gluon loop yields @ 

f dPk^P- ~ gT 2 (2) 
./soft k fa 

because of the enhancement due to the Bose-Einstein statistical factor n B (k). We then find 

n 00 ~ ~ e 2 gT 2 . (3) 
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Of course, we have a cancelation when taking the trace IP^: this is why this diagram is said to be 
vanishing in the HTL framework. The important point is that the vanishing of the HTL contribution 
occurs in this diagram at the order e 2 gT 2 while the calculated soft fermion loop contribution is of 
order e 2 g 3 T 2 . Therefore, non-leading terms (in the HTL sense) in the gluon tadpole may be as large 
as the soft fermion loop contribution. This is indeed the case and these non-leading terms will play 
a very important role. A similar discussion may be given concerning the fermion tadpole Fig. |l|-(c) 
which should however be suppressed by a power of g compared to Fig. [j]-(d) because of the lack of 
Bose-Einstein enhancement at soft momenta. 

In terms of physical processes Fig. ||-(a) involves the first three amplitudes shown in Fig. ||| (and 
their crossing symmetric ones) where the photon is radiated off either by a soft quasi-fermion (a) 
or by a hard on shell fermion (b and c): the imaginary part of Fig. |l]-(a) is in fact constructed 
from squaring Fig. 0-(a) or taking the interference of (a),(b) and (c). On the other hand "squaring" 
the diagram Fig. |||-(b) or (c) yields the soft fermion tadpole diagram of Fig. @-(c). In all cases 
though, the interactions between hard particles in the plasma are mediated by soft fermion exchanges. 
The absence of processes where a gluon is exchanged between the partons in the plasma, which are 
intuitively expected to be relevant, seems to indicate that the calculation based on the soft fermion 
loop is not complete. The expected amplitudes are, in fact, contained in Fig. |l|-(d) and are obtained 
when "cutting" through the effective gluon and the effective 4-point function: these are shown in 
Fig. |2|-(d) and (e). In this set of diagrams the hard gluon line can be replaced by a hard quark line, 
since the radiating quark can be scattered by any kind of parton in the plasma. The physical process 
involved is the bremsstrahlung of the photon by a hard quark which scatters in the plasma via the 
exchange of a soft gluon. The calculation of the diagram of Fig. 0-(d), in thermal field theory, will be 
the main purpose of this paper and it will be found that indeed its contribution to Im IF" M is at least 
of order e 2 g 4 T 3 /qo and that, in the case of quasi-real photons, of virtuality Q 2 /q 2 , <C 1, it is further 
enhanced by powers of 1/g. 

To be complete, one should also mention that the soft fermion loop contribution to real photon 
production has been calculated [||, ^ and the result is as in Eq. (jl]) with the added feature of logarithmic 
mass singularities due to the emission of the photon collinear to the hard massless quark. These 
singularities are regularized by giving the quark its asymptotic thermal mass 

2.2 Semi— classical approach: LPM effect 

The diagrams of Fig. |^-(d) and (e) have been considered in a series of papers by Cleymans, Goloviznin 
and Redlich [[?], ||, |9[. In fact, the purpose of these works was to study the effect of multiple rescattering 
of the radiating quark in the plasma: they find that even when working with massless quarks the rate 
of production of on-shell photons is finite due to rescattering effects, and, more importantly, they 
show that the shape of the photon spectrum is modified, at small energies, and becomes 1 / ^/qo rather 
than the usual logarithmically divergent bremsstrahlung spectrum in l/(/o : this is an illustration of 
the Landau- Pomeranchuk-Migdal effect . The results of @, |9| are obtained in the semi-classical 
approximation (the quark follows a classical path) and it is assumed that the interaction of the quark 
in the plasma is as at T = 0, the only thermal input being the Debye mass, introduced by hand to 
regularize the forward singularity of the quark scattering amplitude. With the LPM effect taken into 
account their result yields ImrP^ ~ e 2 g 2 T 2 (T/qo) 1 ^ 2 , whereas without the LPM effect, the case we 
concentrate upon below, one obtains ImIL" M ~ e 2 g 2 T 3 /qQ. We will come back to a detailed comparison 
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of this result with our calculation after our results have been presented. 



3 Generalities of the calculation 

This section is devoted to explain the notations and to give general results to be used later whatever 
the virtuality of the emitted photon is. 

3.1 Some notations and preliminaries 

The invariant rate of production, per unit time and per unit volume of the plasma, of a real photon is 
given to first order in a {i.e. this formula takes only into account the emission by a quark of a single 



photon, but does not consider the possibility of emitting more than one photon) by [20]: 



R s qo d^k = ~ (2^ ns(9o) ImIF ^- ^ 

For a lepton pair of invariant mass squared Q 2 the rate, integrated over the leptons kinematical 
variables, is 

dN a 1 . . 

It is worth noticing that these formulae, although giving only the dominant term in a of the quark, are 
exact to all orders in a s . We will need to compute the imaginary part of the photon self-energy, which 
is given at lowest order in a s by the two loop diagrams of Fig. ^, to which one should also add the 
diagram with the self-energy correction on the upper fermion line (not represented). These diagrams 
are an explicit representation of the gluon tadpole of Fig. ffl-(d) and, as explained before, they will be 
calculated beyond the HTL approximation. In the following, a summation over the transverse (T), 
longitudinal (L) and gauge dependent (G) modes in the gluon propagator is always assumed. 

The notations for the momenta are explained in Fig. ||. In the kinematical domain we consider, 
the fermion in the loop has a hard momentum, which means that the components of its momentum 
are of order T, and the exchanged gluon as well as the external photon have momentum components 
much smaller than T. 

3.1.1 The "retarded-advanced" formalism 

In order to compute the needed imaginary part, we will proceed by using the "retarded-advanced" 
formalism 22, 23], which is one of the simplest versions of the real-time formalisms. This formalism 



is obtained from the more common "1-2" formalism [24, 26, 27] by a change of basis. The result 
of such a transformation is a formalism with a new 2x2 matrix propagator and new vertices. A 
nice feature of the "retarded- advanced" formalism lies in the fact that it takes into account, in an 
optimal way, the Kubo-Martin-Schwinger identities valid at thermal equilibrium, so that among the 
four components of the matrix propagator only two are non zero; more precisely, the matrix propagator 
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we get for a real scalar field is (for fermions or gauge fields, one has to take into account also the Dirac 
or Lorentz structure): 

A(P) = f A »« P) ° ), (0) 



A„(f>) 

where the retarded and advanced propagators are given by 

Another feature of this matrix propagator is that it no longer contains any information relative to 
the thermal equilibrium. Of course, as said before, the change of basis for the propagator has as a 
consequence a change in the vertices to be used, and the price to pay for these simple propagators 
is somewhat intricate vertices. For example, for a QED-like vertex with entering momenta P, Q and 
—R, Q being the photon momentum, we have: 

g AAA (P,Q,-R) = g RRR (P, Q, -R) = 0, 

g RRA (P, Q, -R) = g ARR (P, Q, -R) = g RAR (P, Q, -R) = e, 

9 R aa( P iQi- R ) = ~ e ( 1 + n B (qo) - n F (r )), 

9 A ra( P iQi- R ) = -e (1 - n F (p ) -n F (r )), (8) 

where n B (qo) = 1/ (^e^ qo — lj, n F (po) = 1/ (^e^ Po + and f3 = l/T is the inverse temperature. 
Therefore, in this formalism, all the thermal information is put into the vertices. The thermal content 
of an n-point function is represented by the external R/A indices. Inside a loop a sum is made over 
the internal thermal indices with appropriate internal propagators joining them. 

3.1.2 Effective soft gluon propagator 

Since the exchanged gluons can be soft in our diagrams, we will need to use effective propagators taking 
into account the Hard Thermal Loop corrections |2|, [|. We choose to use a covariant gauge in which 
the effective gluon propagator is decomposed into its transverse, longitudinal and gauge components: 

-D^(L) = P^(L)A T (L)+P^(L)A L (L) + £P^(L)A G (L), (9) 

where £; is the gauge parameter, and with 

T.L 



A 1 (L) 



A°{L)^j- 2 . (10) 

The appropriate choice of Iq ± ie is implicitly understood when considering R, A propagators. The 
tensors PJjf {L) are the transverse, longitudinal and gauge projectors whose explicit expression are 



25, 29 



P£"(L) =^- — 

1 K z 



P^(L) = U^U U + 

u 



k? L 2 

T {IT V 

= -PS" (11) 
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where 7^ = g^ u — U^U V ', = ■j flu L u , and where U is the mean velocity of the plasma in the current 
frame. As is usual, in the following we will work in the rest-frame of the plasma in order to have 
simpler expressions. 



The explicit form of the functions II Ti (L) are|28|, 29 1 

■•- x(\ — x 2 ) 



U T (l,l ) = 3m 



+ 



n L (Uo 



3m„ 



2 

(1 -x 



In 



x + 1 



4 \x 

„2\ 



1 



, x(l — x z ) , fx + 1 
In . 

1 



(12) 



where x = lo/l and m 2 = g 2 T 2 [N + N{/2]/9 is the gluon thermal mass in a SU(N) gauge theory with 
Nf flavors. Since these functions depend only of x = lo/l, we will often simplify the notation by simply 
indicating x in their list of arguments. 

It will be useful to introduce the spectral functions 



Pt.l(Mo) 



L 2 -U T JL) 



L 2 -Tl T . L (L) 



(13) 



The properties of these spectral functions are closely related to the analytic structure of the gluon 
propagator. In particular, this propagator possesses two thermal mass-shells (a transverse one and 
a longitudinal one) above the light-cone, determined by the equation L 2 = U TL (L), and to which 
correspond delta functions in p TL for L 2 > 0. When L 2 < (i.e. \x\ < 1), the self-energies TL TL (L) 
acquire an imaginary part due to the logarithm (a phenomenon known as Landau damping), so that 
the corresponding contribution in p T L is: 



-2hnH T JL)\ } 



L 2 -ReU TL (L) 



+ (lmU T JL)\ R 



(14) 



For completeness we give here the value of the imaginary parts to be used later: 

3TTm 2 



(15) 



lmU T (L)\ R = —^(1-^) 
ImU L (L)\ R = -2ImU T (L)\ R . 

3.1.3 Resummed hard fermion propagator 

Since the calculation done with massless fermions leads to collinear divergences |i"5]| , it has been 
proposed in [17] to cure these singularities by taking into account the thermal mass even for hard 



particles. Without any approximation, this resummed propagator reads [29, [3(J, 31]: 



S(P) 



7 V - p ■ 7 



+ 



7 + p ■ 7 



(po - S (P))- || p- S(P) || +ie (po - S (P))+ || p- E(P) || -ie. 



(16) 
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where, at the one loop order, 



So(P) 



E(P)=p 



M 2 
4p 

M 2 r 

oo 



2/; 



1 



a; + 1 
x- 1 



x + 1 
1 



X 



(17) 



with x = po/p, and where = g 2 C F T 2 /A stands for the asymptotic fermion thermal mass. 

It is well known that this resummed propagator possesses two mass-shells (usually denoted by + 
and — ), the energy of which can be approximated for hard p by: 



W+(p) « : Jp 2 + M 2 Q 



w- (p) P (l + 2 exp (-4p 2 /^ - 1 



(18) 



Moreover, the residue associated to the — pole is exponentially small whereas the residue associated to 
the + pole is approximately 1. This means that in a momentum integral involving such a propagator 
and where p is hard, we need only to take into account the + pole. In such circumstances, the 
resummed fermion propagator can be simplified to get: 



S(P) 



Pol ~ u + (p)p • 7 
pl - w%(p) +ie 



where 



3.2 Expression of ImFF 



S(P) 



P 2 -M^ + ie' 



(19) 
(20) 



Armed with these tools, one can now evaluate the imaginary part of the diagrams represented on 
Fig. 0. The imaginary part of the photon self-energy is related to the retarded and advanced ones by 



2iImIP, 



IP 



Mr 



IP 



(21) 



In order to simplify the expressions, we will forget for a while all the color and group factors, and 
reintroduce them only at the end of the calculation. 

Let us begin by the self-energy correction of Fig. |3|-(b) (keeping in mind that exists a second such 
diagram) . The contribution of this diagram to the retarded self-energy of the photon is || ^] : 



self 



+ 



d n R 
(2vr) n 

\ ~ n F (r ) 



[S R (P)-S A (P)} [S R (R)] 2 Tr£ R 
S 2 R (R)TTE R -S 2 A (R)TrS A ]S A (P) 



(22) 
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where the notation TrU a with a = R, A stands for: 

TrE a = Tr ( 7 „# [-iZ a (R)] frff>) 



(23) 



A similar relation is easily derived for the advanced self-energy. The one-loop fermion self-energy 
T, a (R) is calculated in a covariant gauge with the effective soft gluon propagator of Eq. (||). The 
thermal distributions can be arranged in a simple way in order to give: 



ImlFV(Q) 



self 



2 2 

; e 9 



d n R 
(27r)™ 



[n F (r ) - n F (p )] 



d n L 

(271-)™ 



A T R L > G (L) - AT' G '(£)] P£JL) Tr; 



self 
p<7 



x [5 fl (P) - S A (P)] [S R (R + L)- S A (R + L)] 



x [n B (l )+n F (r + l )]V 



(24) 



for the transverse (T), longitudinal (L), and gauge (G) contributions corresponding to the propagator 
of Eq. (H), where V denotes the principal value and Tr^ f stands for the Dirac's trace 



Trlf = Tr [fprfa{H + t)l P Hl^] ■ 



(25) 



The reason why we kept only the cut (b) of Fig. |3|-(b) is that when Q 2 > (which is the case of 
interest if one wants to study the production rate of real photons or of dileptons) , the cuts (a) and (c) 
are not allowed by kinematical constraints. Indeed, an on-shell fermion can emit by bremsstrahlung 
only photons with a virtuality Q 2 < 0. 

The same kind of calculations can be carried out for the contribution of the vertex correction of 
Fig. ^-(a) to the retarded self-energy of the photon |2l], §: 



-*n%(Q) 



vertex 2 

R 



d n R 



ppa I j \ rp vertex 
T.L.G\ > P<J 



+ 



2 - n F (Po) 



(2tt)™ t ' l ' g 

(KnTiP Q> -R)s*(P) - v AR f(P Q, -R)s A (P)) s R (R) 

(V AR L f(PQ,-R)S R (R) - V^f(P,Q,-R)S A (R)) S A (P)} 



(26) 



again corresponding to the transverse (T), longitudinal (L), and gauge (G) terms of Eq. (0), where 

r-Q, vertex 



, ,, T contains the Dirac's algebra factors: 

Trr tex _ Tr [ 7 ^ 7p( ^ + fy^tf + ^) 7(j ^ 



(27) 



and the functions V^' (P, Q, —R) are the scalar part of the one loop q — q — 7 vertex, where the 
gluon line of the loop is made of an effective propagator. These are equal to 



V^°(P,Q-R) 



d n L 
(27t)™ 

A T R L ' G (L) - A T A L - G (L)} S a (P + L)S^(R + L) 
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+(\ - n F (ro + *o)) [S R (R + L)-S A (R + L)] S (P + L)A^ G (L) 



1 



+(5 " u f(po + k)) [S R (P + L) - S A (P + L)\ S P (R + L)A. 



7 

T,L,G 



a 



(28) 



Inserting this expression into Eq. ( |26[) , and calculating in the same way the advanced self-energy, one 
obtains the contribution of Fig. ||-(a) to the imaginary part of the photon self-energy: 



1 2 2 
-e q 

2 y 



d n R 

(27T)™ 



[n F (r ) -n F (p )] 



d n L 

(2TT) n 



A T R L - G (L) - A T A L ' G (L)} P^ G {L) Tr™*- 



x I [S R (P] - S A (P)] [S R (R + L) — S A {R + L)] 
x [n B (l )+n F (r + l )}V 



1 



R 2 -Ml 



V 



1 



+ [S R (R) - S A (R)} [S R (P + L)-S A (P + L)} 
x [n B (l ) +n F { Po + 1 )}V ' 



M 2 

00 



P 2 -Ml 



V 



1 



(it! + L)2 - Ml 



(29) 



For the same kinematical reasons as in the self-energy diagram, the cuts (a) and (c) do not contribute 
to the production rate of real photons or lepton pairs. We notice at this point a high symmetry between 
the self-energy correction and the vertex correction. Two kind of terms appear in the latter and can 
be interpreted in terms of cut diagrams: the first term in the curly brackets above corresponds to cut 
(6) and is to be combined with Eq. ( |2"4| ) while the second term should be combined with the other 
self-energy correction mentioned above but not made explicit. Both classes of terms give identical 
contributions to the photon production rate, so that in the following we consider only the first term 
of Eq. (^9|) and Eq. (^4|), and take the other terms into account simply by multiplying by an overall 
factor 2. 



3.3 Reduction of the traces 

In order to make some simplifications more obvious, it is useful to compute the Dirac's traces Tr^ f 
and Tr™ rtcx in such a way that the invariants P 2 , R 2 , (P + L) 2 and (R + L) 2 appear whenever possible. 
We get: 



r self 



AR 2 Q p R a - AQ 2 R p R a - g pa (r 2 {R 2 - Q 2 ) + 2R 2 Q ■ L - 2Q 2 R • L 



and 



(30) 



rp ^ vertex a 



2(R 2 + (R + L) 2 )P p Q a - 2{P 2 + {P + L) 2 )R p Q a 
+2L 2 (R p R a + P p P a ) - AQ 2 R p P a 
+9 P a (P 2 R 2 + {P + Lf{R + Lf 

-L 2 (P 2 + R 2 + (P + L) 2 + {R + L) 2 - Q 2 - L 2 ) 



(31) 
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At this point, it is worth noting that we dropped terms proportional to L p or L a , since they will 
disappear in the contraction with the projectors P£ a L - 

L p P£° L (L)=0. (32) 

Therefore, the expressions obtained above for the Dirac's traces should not be used to compute the 
gauge dependence of the production rate, since L p PP a (L) / 0. In order to verify the independence 
of this rate with respect to the gauge parameter £, it is simpler to go back to the original expressions 
Eqs. (p4j), (p9|) and show a compensation of the gauge dependence between the self-energy and the 
vertex contributions, exactly in the same manner as at T = 0. 



3.4 Kinematics 

We now take into account the constraints provided by the discontinuities: 

S R (P) - S A (P) = 2ire(p )5(P 2 - A&) 

S R (R + L)-S A {R + L) = 2^e(r + l )5((R + L) 2 - M^); 



(33) 



this will enable us to perform two of the momentum integrals "for free". Moreover, as already men- 
tioned, the components of R and P are considered to be hard, whereas the components of Q and L 
are considered to be much smaller than the temperature T. Therefore, po and ro + Zo are of the same 
sign: 

e(po)e(r + Zo) = I- (34) 
From the constraint 5(P 2 — M^), we extract the value of tq: 



q ±Jr 2 - 2qr cos 9 + q 2 + 



a 2 M 2 
<cos# + !-(l-cos 2 0) + -iH 

2r 2r 



(35) 



where 9 is the angle between the spatial components of R and Q. Since the possible two signs for 
ro lead at the end to equivalent contributions, we will consider only the positive solution for ro, and 
multiply when necessary the result by a factor 2 to take into account the other solution. From the 
constraint 5{{R + L) 2 — M^,), we obtain the value of cos#', where 9' is the angle between the spatial 
components of R and L: 



cos 9' 



(r + l ) 2 -(r 2 + l 2 +Ml) 



2rl 

qo — q cos 9 
I ' 



(36) 



Moreover, we must impose cos#' to be in [—1,1], which can give additional constraints between 
Zo, Z, qo, q and cos 9. When needed, the variable cos 9", where 9" is the angle between Q and L, will be 
obtained by: 

cos 9" = cos 9 cos 9' + sin 9 sin 9' cos (f>, (37) 
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where (f> is the angle between the projections of Q and L on the plane defining the spherical coordinates 
basis. Given these relations, the angular variables we choose as independent ones are 6 and <p, 0' being 
extracted from the constraint Eq. ( |36| ) and 6" being given by Eq. (37). 

With such a choice of independent variables, we can calculate approximate values for the various 
denominators entering the rate: 



R 2 - M^ = 2Q ■ R-Q 2 « 2q r 



1 



/If* q: 



(38) 



2- 



2tt 



Q 2 - 2Q • R - 2Q ■ L 



2tt 



2q r 



Ml , Q 2 



1 1/2 ' 



(39) 



(the reason why we performed the integration over <p at this point is that at the dominant order 
(P + L) 2 — M^, is the only quantity in which appears). 

It is worth noticing that with = and Q 2 = 0, the denominator R 2 — exhibits a collinear 
singularity when 9 = 0, i.e. when the photon is emitted collinearly to the hard fermion. Looking at 
the previous expressions, we see that it is natural to introduce the new variables: 



u = 1 — cos I 



M. 



off 



2 Q 2 r 2 
M + — 

% 



(40) 



It is now clear that the collinear divergence is regulated by the effective fermion mass M c g which acts 
as a cut-off in the integral over u. Therefore, the order of the result depends on the order of this 
cut-off, a question which will be discussed in the following paragraph. 



3.5 Enhancement mechanism 

A priori, we are faced with four kinds of angular integrals over the variable u: 

2 2 



h = 



du 



R 2 -Ml> 



h = 



du 



du 



(R 2 -Miy 



h 



du 



(R2-Ml)((P + L) 2 -Ml 



(41) 



Since the order of R 2 — or (P + L) 2 — is q$r, the orders we expect naively for these integrals 
are: 



O 



\qor 



o 



1 



(42) 
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and this would be right, up to logarithmic factors, if the collinear singularities where really logarithmic 
in the four integrals. This is actually the case for I\ and I2 but obviously not for ^3, which exhibits a 
power- like collinear singularity. Looking at the form of R 2 — M^, we have for the order of I3: 



where u* ~ M 2 s /r 2 . Therefore, the actual order of I3 depends on the order of the physical cut-off in 
the u integral, this dependence being a power-like one instead of a logarithmic one. Moreover, this 
integral is enhanced with respect to its naive order when w*<Cl, which occurs when the virtuality of 
the emitted photon is small. 

As for I4, since the hard momenta R and P + L only differ by soft quantities, the denominators 
R 2 — M^j and (P + L) 2 — are zero almost simultaneously. Let us quantify more precisely this 
assertion; to explain roughly what happens in this case, it is sufficient to write R 2 — ~ qor(u + u*) 
and (P + L) 2 — ~ qor(u + u*'). Again, we have u* ~ M 2 s /r 2 . Moreover, it*' — u* is independent of 
the value of M 2 s /r 2 and is of order L 2 /r 2 <C 1. Therefore, we get the following orders for the integral 



Like I3, this integral is enhanced when «* « 1, The reason why I4 behaves much like I3 whereas I4 
contains two simple poles instead of a double one in the case of I3 lies, of course, in the very close 
vicinity of these two poles. In fact, when u* is smaller than the distance between the two poles, the 
order of the result is controlled by this separation which is of order L 2 /r 2 . On the contrary, when the 
regulator u* is larger than this separation, everything happens as if we had a double pole. 

Recalling now that the fermionic thermal mass is of order gT, we should distinguish various 
cases according to the value of the virtuality Q 2 (see Eq. (|40|)): 

(i) Q 2 jq\ <C g 2 : in this case, the virtuality of the emitted photon does not play any role, and can 
be completely neglected, since M c g « M^. This case will be studied first since it is the simplest 
one, and will provide a basis for the case where the virtuality can no longer be neglected. In 
this case, the integrals I3 and I4 are enhanced whereas I\ and I2 have their "normal" order, 
therefore we can neglect every occurrence of I\ or 1^- 

(ii) g 2 < Q 2 /ql < 1 : in this case, the virtuality must be taken into account since its effect is at 
least as important as the effect of the fermion thermal mass. We still have u* < 1, so that the 
remark made in (i) on the enhancement of ^3 and I4 remains true. This section will make an 
intensive use of the results obtained in the case (i). 

(iii) Q 2 /qQ ~ 1 : here also, the virtuality of the photon plays an important role. Moreover, this case 
is very different from the previous two since none of the Ij's is enhanced, and therefore all must 
be taken into account. Another reason why this case is very different lies in the fact that there 
is no more a hierarchy between the various powers of Q 2 /q 2 ,; as a consequence, all the powers of 
this quantity should be kept in the expansion of the numerator. 




(43) 




(44) 
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It is interesting to compare the status of mass singularities in our calculation with results obtained 
in previous works. Collinear singularities and their cancelation have been studied in (non-resummed) 
perturbation theory at two-loop order in QED/QCD J32| , |33| ] and three-loop order in the A</> 3 model 
[ 34 1 : in all cases the problem considered was the decay rate of a heavy particle into massless particles 
and a complete cancelation of collinear singularities associated to the massless particle propagators was 
found when summing over real and virtual diagrams and over all thermal processes. The equivalent 
problem here is the case Q 2 7^ but = and indeed our expressions are regulated by the photon 
virtuality: in the R/A formalism there is no need to distinguish between "real" and "virtual" diagrams 
as they are all included in rather compact expressions (see e.g. Eq. ( f45[) below). It is the masslessness 
of the "decaying" particle (Q 2 = 0) which generates the collinear singularity of interest here as it 
would in the previous studies if we had let Q 2 — > 0: indeed the two-loop correction to the invariant 
rate R was found to be ~ e 2 g 2 T 2 ln(T 2 /Q 2 ) [35], which is indeed logarithmically divergent when the 
photon virtuality vanishes. This softer singularity (compared to the "power- like" divergence above) 
can be understood because no HTL resummation was performed on the gluon propagator (L 2 = 0) 
and therefore no Landau damping was included. 



4 Production of quasi real photons: Q 2 /q$ <C g 2 

Let us now specialize to the production rate of almost real photons, which corresponds to the case 
denoted by (i) at the end of the previous section. From a technical point of view, this case is the 
simplest one. Nevertheless, it will serve as a basis for the production rate of photons having a small 
virtuality, since this latter case will appear as a generalization of the present case. 

4.1 Expression of ImlP ^ 

By using the fact that Q 2 = and the fact that P, R are hard and Q,L are much smaller than T, 
it is possible to greatly simplify the expression of the imaginary part of the photon self-energy. In 
particular, we get some partial cancelations between the contributions of the self-energy and of the 
vertex. The remaining dominant term is: 

/d n R ( d n L 
( 27r )n-l J ( 2 ^)n-l g° W F(P0) n B (h)p TL {l,h) 

xe( Po )e(r + l )8(P 2 - M^S [(R + L) 2 - M 2 ^ 

L 2 

x (R p R a + P p P a )P£ L {R2 _ M^yjp + L)2 _ Ml) . (45) 

where we took into account the fact that the integral I4 is enhanced. We have also simplified the 
statistical factor n F (ro) — n F (po) to qo n ' F (po) which is justified within our kinematical approximation. 
The fact that the first non vanishing term in the numerator is proportional to L 2 <C T 2 shows clearly 
that ImIT M ^(Q) is zero at the HTL order. 
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Using Eqs. (p5|), (p6|), the expansion of the terms appearing in the numerator gives 

2 

L 2 (R p R a + P p P a )P?° 



7 1 < iM 



L 2 {R p R a + P p P a )PP° 



(46) 



Using now the constraint — 1 < cos #' < 1 (Eq. (|3q)), we see that the gluon momentum L is constrained 
to satisfy Iq/1 € [—1, 1], i.e. L 2 < 0. This means that the production rate of real photons is dominated 
by the exchange of a space-like gluon (i.e., by the Landau damping part of the gluon spectral density). 
Eq. (f45|) can be then be put into the form 



ImlT%(Q) 



e 2 g 2 1 
T 8vr4 go 

2 



2\2 



drn' F (r) j l*dl / dx n B (lx) p T L (l,lx)(l — x z ) 



x 



o 



"27^" 



u + 



du 



2r 2 



+ 



"2T 2 " 



+ 



1/2 ' 



(47) 



where we denoted x = and where the symbol (— 1) T denotes an extra minus sign in the transverse 
contribution. We have introduced some cut-offs r* and I* at a scale intermediate between gT and T, 
since we a priori restricted r to be hard and I to be negligeable in front of T. We will see later that 
these cut-offs can be taken respectively to and +oo, without modifying significantly the result. 

An alternate expression for Im EP ^{Q), which will be needed later, can be obtained in the following 
way. Using the 5— function constraints one easily derives 

1 1 1/1 1 



R 2 -Ml{P + L) 2 -Ml 2Q-L\R 2 -Ml (P + L) 2 - M 2 



+ 



(48) 



Then, noticing that the change of variables P + L — » —R leaves the integrand in Eq. (|45|) invariant 
it is legitimate to replace in this equation [(R 2 - M^)((P + L) 2 - M^)]- 1 by [-Q ■ L(R 2 - M^)}- 1 . 
We are then lead to the expression [] 



Imn%(Q) 



eVl 

c+1 



r 2 dm' (r) 



2\2 



x I dl dx n B (lx) p T L (I, lx)(l — x ) 



i 



du' 



VT^T! 7 4M^ + £2(1 _ x 2y ' 



(50) 



a An alternative method to obtain the same result directly from Eq. (|47[) is to perform a few changes of variables 



lM a 



v 

1 /Qo 

2 V a 



a; 2 t + 



Q 

do 



l 2 (l-x 2 ) 
2r 2 

with «o 



Ml 
2r 2 



l - ^T^" 



(49) 
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where the relation between the u' and u integration variables is v! = 



-8r 2 u/L 2 . 



4.2 Reduction of the expression 



The expression of Eq. ( |50| ) can be simplified by performing the angular integration over the variable 
v! . This integral is elementary and gives: 

dv! 1 _ 2 tanh-y-L7(4M^ - ^ 



v T^?4M2 3 +Z 2 (l-x 2 )n' ^-L 2 (4M^ - L' 2 ) 

At this point, we notice that the integral over r can now be factorized in order to obtain: 

+00 +00 

drr 2 n' F {r) ~ / drr 2 n' p (r) = — . (52) 

r* 

The reason why we can set the cut-off r* to zero lies in the fact that the function to be integrated is 
peaked around r ~ T, whereas r* <C T. Indeed: 



7 

/ 



drr 2 n' F (r)^-^L. (53) 



Moreover, in order to further simplify the result, we will approximate the Bose-Einstein weight by 
n B (lx) ~ T/lx since the momentum / is assumed to be much less than T. 

Finally, if we recall that the spectral density to be used here is given by Eq. ( |l4| ) since L 2 < 0, and 
if we introduce some convenient dimensionless quantities: 



-L 2 



ImIT 



I T ,l( X ) - y2 
lv± oo 

~ Re II™ T (x) ., 

R T ,A X )= M2 ■ ( 54 ) 



the imaginary part of the real photon self-energy can be written as: 

„ , s , s e 2 g 2 N n C T? T 3 } dx ~ , s 

07T Z qo J X 



dw ^/w/jw + 4)tanh~ 1 y /w/{w + 4) 
(ti; + fl r , i (x))a + (J, J> (x))2 ' 



where we have re-introduced the required color and group factors. 
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4.3 Reduction using sum rules 

Instead of performing first the angular integral over u to reduce Eq. ([17]) to the two-dimensional result 



of Eq. (pq), it is possible to use sum rules to obtain the result as a one-dimensional integral. 

Let us first recall that sum rules come in this context from the spectral representation of the 
resummed propagator, which can be written for the scalar-like part of the gluon propagator as: 



+00 



ll -I 2 -U TL (l,l ) +ie 



E 



dE 2^ p t,l{ 1 ^ E ) Y 



E 2 + 



(56) 



By taking the imaginary part of this equation, keeping in mind that p T L are real functions, we can 
then obtain: 



dx 
2^ 



x p TL (l,lx) V 



1 



l 2 (y 2 -l)-ReU T Jy) 



P {y 2 _ 1} _ Reu T ,(y)) + fimn T , (J,) 



2 • 



(57) 



where we denoted Iq = ly and E = Ix. By taking now special values of y, one easily obtains: 

-I- no 

dx p T (l,lx) 



With y = 



-00 
+00 



-00 
+00 



With y = 00 



With \y\ > 1 



2-7T X 

dx p L (l,lx) 
dx 



1 



1 



x 



I 2 + 3m 2 



-00 
+00 



2tt 

dx 
2^ 



x p TL (l,lx) = 
x p TL (l,lx) V 



1 

l 2 



l 2 (y 2 -l)-ReU T AyY 



(58) 



These are the basic sum rules we will use to perform the integral over x. However, the integrals we need 
to perform are on x £ [—1, +1] whereas the above sum rules give the result for an integration over the 
whole real axis. Therefore, it will be necessary to subtract the contributions of x G] — 00, — 1] U [1, +oo[. 
This subtraction is easy to perform since in that region, the spectral functions p TL (l,lx) are made of 
delta functions corresponding to the thermal mass-shells of the gluon. More precisely, we have: 



/ 



]-oo,-l]U[l,+oo[ 



T.L 



f 



I 



f 



I 



(59) 



where Z TL stands for the residue of the spectral function on the mass-shell, whose energy is u> TL (l). 
Explicit forms of these residues can be found in Ref.Q; for later use we record here the asymptotic 
limits 



3 

Jt{l -» 00) ~ I 2 + -m\, 



uj 2 (1 -»■ 00) ~ V 



1 + 4 exp 



2l 2 
3m 2 



16 



Z T (l — > oo) 



3m 2 
U 2 



In 



4Z 2 
3to 2 



Z L (l -> oo) 



2Z 2 
3m 2 



4Z 2 

— K exp 



2Z 2 
3m 2 



Z T (Z -+0) ~ 1 



5 
Z 2 



u, 2 (Z^0)^ + ^ 2 
2Z 2 



5mj 



ZJZ^0)~1 + 



5m 2 



(60) 



As well, we can introduce in a natural way a magnetic mass for the transverse gluon by assuming that 
it modifies the static behavior in such a way that: 



n T (o) = m 2 n 



(61) 



which, thanks to Eq. (p7|), is equivalent to: 

J £-xp T (l,lx) = p, 



-oo 
+oo 



fix /9 T (Z, Zx) 
2ir x 



I 2 + m: 



2 

mas 



(62) 



These sum rules may now be used to reduce Eq. (|5Cj) down to a one-dimensional integral. To do 
this, we first cast Eq. (|50|) into the form 



2„2 rp3 



eg 



Imn%(Q) «(-l) T .., , 

127T 2 go 



1 



i* +i l 
Idl \ dx p T L (I, Ix) 



du' 



o 



-i 







_±_|_ (x 2_ 1)+2x(y 2_ 1) + (y 2_ 1)2 



i i 

+ 



x + y x — y 



(63) 



where y = y/l + 4M^/I 2 m'. The application of sum rules to evaluate Eq. (|6^) is now completely 
straightforward. We consider each of the terms in curly brackets in turn, and examine first of all the 
transverse contribution. 

For the first term, we have a contribution due to 



I* +1 l 
Ki = -2 J Idl J dx p TL (I, Ix) J — 



du' 



1 2 



-i 



d! 



v/Z 2 + 4M 2 



In 



u'y/1 - u' 2y 2 x 
o 

(Z + ^ + 4M£) 2 ~ 



x 2 -l) 



4M 2 



(1 - Z T ) - V 



Z. 



P + m*^ u 2 



(64) 



where we have introduced a magnetic mass m n 



g T via Eq. (|62|). Thanks to Eq. (pQ), this integral 



is finite for Z — » 0, and in fact for M m ~ (7T we can drop m 



mag 



fT, but not if ~ g T or below. 
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In the following we use the HTL result, 
magnetic mass. It is also easy to see from Eq. 



gT, so that we can drop the dependence on the 
that since 



m„ 



dl 



I 



m a \ i 



I* 



(65) 



we can take I* —* oo by introducing a negligeable contribution. 

For the remaining terms of Eq. (|63|) , we find the sum rules lead to 



i* +i l 
K 2 = -2 f Idl [ dx p (l,lx) I 



du' 



1 



o -1 

I* l 
dl { du 



o 



AMi 



2x(y 2 - 1) + (y 2 - If 



Ren T (y) 



1 1 

+ 



x + y x — y 



I J VT^ 7 4Ml + l 2 u' \ T AM^-u'^l -I 2 ) AMI - u'ReU T {y) 



(66) 



We note that no magnetic mass term arises for this contribution. In this form, using Eqs. (fjB), it is 
straightforward to verify that Eq. (|66| ) is finite as I —* 0, as there is a cancelation between the two 
terms in curly brackets. As well, as was the case of K\, we can take the limit of the cutoff /* — > oo by 
introducing a negligeable contribution. To carry out the v! integration in Eq. (|66|) it is convenient to 
add and subtract a term representing the small I behavior: 



oo 1 

dl 



du' 



AMI 



I J VT^T! 7 AMI + l 2 u' 





m„ 



' AM 2 , - u'{uj 2 - I 2 ) AM 2 , - m 2 u' 



Ren T (y) 



m„ 



AM 2 , - u'Re II T (y) AM 2 , - m 2 u' 



■ (67) 



The integration over u' can now be carried out, but the explicit results depend on the signs of 1 — 
4M^/Ren Ti (x), 1 - AM1 /(uj 2 tl - I 2 ) and 1 - AM^/m 2 , where \ = y/l + Mi^/F. The quanti- 
ty 1 - 4M^/Ren Ti (x) is plotted on Fig. ^, for both the transverse and longitudinal case. Since 



1' 



Re II T L (x 



/I), and since x is varying from 1 to +oo, these curves give us all 



the signs we need if we recall also the limits lim IE ^ 1 + 1 — AM^,/ReIl T (x) = 1 — 8M^,/3m 2 and 
lim a ._^. 1 + 1 — 4M^/ReII i (x) = — oo. The case 3/8 < M^/m 2 is the most interesting one physically 
(e.g., an SU(3) gauge theory with less than 10 light flavors), and fortunately the simplest one. We 
find for the total transverse contribution the result: 



lmWJQ)\ 



e 2 g 2 T 3 



oo 

dZ J K Itanh- 1 - 
3vr q J I I I 2 + m 2 x X 



+ 



AM 2 

oo 



1 



l 2 +AM 2 r~AMi 



tan 



Ren T ( x ) 



4M2 



WRen T (x) V 
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Likewise for the longitudinal case we obtain: 



Imn%(Q)| 



eV T 3 Tdl 



>~>K 



3mi 



3ir go J I I W 2 + m ? ^ 2 + 3ml I Y 
o 



tanh 



-i 



X 



+ 



AM 2 



oj 2 -l 2 



oj 2 



tan 



Ren L ( x ) 



4M 2 

oo 



tan 



4M2 



\V Ren iM / 



4M 2 

^ — i 2 



7 



+ 



4M2 



1 



/ 2 + AM 2 l 2 +m 2 W 



tan 



4M 2 i 

oo 



(69) 



M 2 

oo 



There are two other cases that need to be considered in general: 1/4 < M^/m 2 < 3/8 and 
/m 2 < 1/4. These two cases are more intricate because now the required signs depend on the 
value of I. For this reason, the corresponding expressions will not be written here. There is one 
limit of interest, — ► 0, which is included in the latter regime, but this limit is more conveniently 
handled at an earlier stage of the calculation (see Eq. (|55|)). One should note that this splitting in 
three ranges is only a consequence of the use of sum rules and has no physical meaning, since it is 
obvious from Eq. ( |50D that the result is a continuous function of the ratio /m 2 . The reason for this 
can be seen from Eq. (|50|): one must extend the x integration range from [— 1,+1] to the whole real 
axis in order to use sum rules and then subtract the contribution of the region ] — oo, —1] U [+1, +oo[ 



by Eq. (||). While x 6 [-1, +1], the denominator 4M^ + r(l 



)u' remains strictly positive; but 



if x g] — oo, — 1] U [+l,+oo[, one can be faced with a pole in the v! variable, which is dealt with a 
principal part prescription. A by-product of this pole is that the expression that should be integrated 
over I depends on which side of the pole we are. 
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4.4 Discussion of the result and asymptotic behavior 

In this paragraph, we study the behavior of ImIP M (Q) in some limiting cases for and m 2 . We 
work with the compact form Eq. (|5^). In fact, the quantity of interest is 

j TL = *i i tl{x ) [ dw VWWTT)t^VWWTT) 

J x T ^ > J (w + R T Jx)) 2 + (I T Jx)) 2 

which is dimensionless and is therefore a function of the dimensionless quantities Z*/M 00 and mjM^. 
The independence on I*, as well as the infra-red finiteness, can be shown to hold using this formula, 
but since we already did this using the sum-rule formalism we will not repeat the discussion here. We 
thus assume J TL to be only a function of M^jm^ and study, both analytically and numerically, the 
limit of the above expression when the ratio M^jm^ — > in order to better understand the origin 
of the singularities. Of course, the precise value of this ratio is fixed by the number of colors and 
flavors of the studied theory, and is of order 1. Despite that fact, this limit enables one to verify 
what is the effect of switching the fermion thermal mass off, and gives information on the nature of 
the singularities which are regularized by this mass. Some details of the analysis are given in the 
appendix. We find: 

J L ~ In 



Jt ~ Wtt 5 ")- (71) 



These results are interpreted as follows: the thermal fermion mass regularizes the collinear divergence, 
and therefore this kind of divergence, common to both the longitudinal and the transverse cases, 
contributes as one power of ln(m g / M^) . Additionally, in the transverse contribution, there is an 
extra power of that logarithm, which is a consequence of the potential infrared divergence due to 
the absence of a thermal mass for static transverse gluons. This means also that this divergence is 
unexpectedly regularized by the fermion thermal mass. In fact, we can see already these features in 
Eq. flsT]), where the tanh -1 function corresponds to the collinear logarithmic factor and where the 



combination y/— ^{AM 2 ^ — L 2 ) (instead of —L 2 if is not taken into account) in the denominator 
is responsible for the regularization of potential infrared problems. We can also observe in that limit 
the essential difference between a true double pole and the structure of the integral we denoted by I4 
in the section |3.5| (indeed, as seen before, I4 is the integral of interest in our problem, because it is 
enhanced, and because it has the dominant part of the numerator associated with it). If we really had 
a double pole, like in I3, we would have got a factor m^/M^ instead of a logarithm. This point will be 
discussed further when we compare our results with that of |7|, |8], |9| . We have checked that numerical 
calculations based on Eq. ( j70|) are in good agreement with the analytical results in the asymptotic 
region. 

In fact, if the fermion thermal mass is sufficiently small, there can be a competition between the 
fermion mass and an hypothetical magnetic mass, expected at the order m mag ~ g 2 T. If we introduce 
this mass by the transformation of Eq. (|6ll) , then we obtain: 

(i) If m mag < Moo < rn s : 
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(ii) If Moo < rn mag <C m g : 



This competition between the magnetic mass and the fermion mass has been studied numerically, 
which gives the curves of Fig. (||), where are plotted the value of J T as a function of log(m mag /7?i g ), 
for various fixed values of the ratio M^jm^. We observe that the effect of the magnetic mass becomes 
non negligeable as far as m mag > 0.1 -Moo. When the magnetic mass is small enough, the flattening of 
the curves denotes the independence of the result on m mag . 

5 Comparison with the semi-classical results 

The purpose of this section is to compare the results provided by the thermal field theory techniques 
with those obtained via semi-classical methods. 

5.1 Structure of the emission rate 

In the semi-classical approach, one finds 0, |, g, H that the emission rate of soft photons in a 
scattering process can be factorized as the product of two terms: the intensity of emission which is the 
square of an electromagnetic current, and the cross section for the same scattering process without 
photon emission. Our purpose is to make a precise connection between such a result and our result 
Eq. (HI). 

Let us begin by defining the amplitude corresponding to the scattering process without any photon 
emission, i.e. corresponding to the diagrams of Fig. ||-(d),(e) where the photon line is suppressed 
(written here for the scattering of a quark on another quark, but it is very simple to write it also for 
the scattering on a gluon; moreover, like in Eq. (E5j) , we don't write the color and group factors): 



iM(Pi,P 2 ;P{,Pi;L) = -u(P[)(-i 9 r)u(Pi) u(P! 2 ){-igY)u{P 2 ) 

x E P^(L)*A a (L), (74) 

a=T,L 

where Pi, P 2 are the incoming momenta, P[, P'^ are the outgoing momenta, and L is the transferred 
momenta. Of course, only three of these five momenta are independent, and appropriate Dirac delta 
functions will be needed to enforce the momentum conservation. The amplitude squared is simply 
given by: 

|A4| 2 = /lH7 M M7>i)^(7^ 2 7^ 2 ) E P^(L)P b p AL)*A a (L)[*A\L)]\ (75) 

a,b=T,L 



We are now going to rewrite Eq. 45 into a form exhibiting the decomposition mentioned above. In 



order to make the connection with Eq. 45, we first notice: 

A(R p R a + P p P a ) « Tr( 7 "(# + t)rP), (76) 

and 

4e^ ^ 2 fP-e (R + L).e \ 2 

- Ml)((P + Lf - Ml) ~ & p ^ £ \P-Q (R + L)-Qj ' {U) 
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where the sum runs over the polarizations of the emitted photon. We recognize the standard electro- 
magnetic current which appears in such a soft emission |3(| . It is worth recalling that we used the fact 
that the relevant values of u are of order ti*<l. We need also the hard thermal loop contribution to 
the imaginary part of the polarization tensor of the gluon: 



ImD&jL) 



9 



d 4 K 



[n F (k + Z Q ) - n F (ko)] x 

2ttS(K 2 ) 2tt5((K + L f) Tr^fC-ftfC + £)), 



(78) 



where we considered only the quark loop contribution. Of course, the gluon loop would be related 
to the scattering of the quark on a gluon, which has not been written. Moreover, like in Eq. ( |45| ) 
and Eq. ([T5|), we have not written the color and group factors. Concerning the statistical weights of 
Eq. a straightforward calculation gives (the n B (qo) factor should be added in order to reconstruct 
the emission rate defined in Eq. (|j), and the other statistical factors have been taken at an earlier 
stage of the calculation, for example in Eqs. (p4|) and (|29|), before any approximation): 



(79) 



n B (?o) [n B (h) + n F (r + h)] [n F (r ) - n F (p Q )] [n F (fc + l )-n F (k )] 

= n F (r + Iq) [1 - n F (p )] n F (k ) [1 - n F (k + Z )] 



Imn gluon( L ) 



a,b=T,L 



We need also the following identity: 

£ P^(L)P^(L)*A a (L)[*A 6 (L) 

2 Imn c (L)P7(L)P;,(L)^(L)*A fl (L) [*A b (L) 

i,b,c=T,L 

£ Imn a (L)rA«(L)| 2 P; CT (L) 



a=T.L 



- £ Pa{L)P*,(L), 



(80) 



a=T,L 



the last equality being valid in the region where L 2 < 0. 

Given all these relations, it is straightforward to rewrite our result for the emission rate (see Eqs. (0) 
and (|45|)) per unit time and unit volume as: 



dN 



d 3 q 
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(81) 
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where \A4\ 2 is given in Eq. (ffq) . This expression shows clearly that thermal field theory calculation 
leads also to a separation in two factors: the amplitude squared of the similar scattering process 
without photon emission, and a factor which is nothing but the square of an electromagnetic current 
and is specific to the photon emission. These two factors are then integrated over the momenta of 
the unobserved particles, with appropriate statistical weights. Because of this high similarity with the 
semi-classical expressions, this structure can be seen as another evidence to say that the enhanced 
terms we have exhibited are actually the relevant contributions to photon production by a plasma. 

This decomposition is valid only when the photon momentum Q is negligeable in front of the quark 
momentum, so that the amplitude that appears in Eq. ( |8l[) is almost unperturbed by the momentum 
Q. It is worth saying that we used several times the fact that u* <C 1 to give that structure to Eq. (fi5|). 
It means that this kind of factorization will not be valid in the case where the virtuality is important 
Q 2 /lo ~ 1- Moreover, compared to the standard semi-classical treatments J?], |8|, ||, we see that this 
structure remains valid even if we take into account the transverse exchange contribution, which has 
been shown in the previous section to be of the same order of magnitude as the longitudinal one. 



5.2 Factorization limit 

Nevertheless, despite that nice structure, the phase space over which the integration has to be per- 
formed is common to the two factors. Indeed, these two factors share some variables, like the momen- 
tum L exchanged during the scattering. Therefore, from a technical point of view, the factorization 
is not really complete. 

The integral over the angles 9 and of the photon line affect only the factor related to the photon 
emission, and is easily obtained via Eqs. (77), (|38|), ( |39| ) and (|5l|): 



2tt 2 



^ (Pi-e P{-£\ 16vr / —L 2 , / —L 2 

o o ""J, live " pFq) a igyiMi^ tmh ~ V «^ (82) 

As already said, the tanh -1 factor is a remnant of the collinear singularity (when M m — ► 0, this factor 
diverges logarithmically) . 

Then, we see that this emission factor still depends on the variable L 2 , so that the apparent 



factorization provided by Eq. 81 is not a true factorization. We can look at the two limiting cases: 

(i) If L < Moo, we have: 

27T 2 n „ 

*£^-*q) (83) 

Therefore, the coupling between L and M^ in the emission factor will be useful at small L 
to regularize potentially dangerous terms in the infrared sector. This is the reason why the 
transverse gluon exchange is infrared finite despite the lack of thermal mass for static transverse 
modes. 

(ii) If L 3> Moo, we have: 

2tt 2 







Pi-e 
Pi-Q 




pol. £ 
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This limit is known as the "factorization limit" |[7|]. Indeed, in that limit, the L dependence of 
the emission factor is only logarithmic and is therefore negligeable, so that the only relevant L 
dependence is in the scattering amplitude. 

This limit is valid for the longitudinal contribution when m g ^> M^, because in that case the 
relevant infrared regulator is m g (which is the thermal mass of static longitudinal modes) and the 
relevant values of L are also of order L ~ m g 3> M^. In that limit, the longitudinal contribution 
becomes proportional to \n(m^/M^ ). 

Almost the same interpretation can be given for the transverse contribution: we get one factor 
In (m^/M^) from the tanh -1 when m g M^, whereas the second power of that logarithm comes 
from the singular infrared behavior of the scattering amplitude. In that case, the "factorization 
limit" cannot be used for all the values of L since values of L of order are explored due to 
the lack of transverse static thermal mass. 

In conclusion, we can say that the expression Eq. (|8l|), despite its high similarity with semi-classical 
ones, should be interpreted with care concerning the possibility of "factorization". Moreover, this lack 
of factorization at small L is closely related to the possibility of getting a finite transverse contribution. 



5.3 Order of magnitude 

After having seen that our expression for the photon production rate, derived from thermal field theory, 
is similar to the expressions usually encountered in semi-classical methods, it remains to compare the 
order of magnitude of the rate given by the two methods since the calculations are performed in very 
different ways. 

From Eq. (||) and the discussion in the previous section, the differential rate for real photon 
production, integrated over angular variables, takes the form 

iN -1)^, m ,85) 




where one recognizes the typical 1/qo bremsstrahlung spectrum. It is instructive to compare our result 
with that in ^, ||] in the case the LPM effect is turned off. It is found there, neglecting numerical 
constants, that in the semi-classical approximation 

r6 f rpl \ 

(86) 

The structure of this equation can be understood as follows: the T 2 /M^ piece is arising from the 
collinear radiation of the photon by the hard quark, while the factor ln(T 2 /m 2 ) reflects the divergence 
of the quark scattering cross section in the plasma, screened by the Debye mass m g . Eqs. ( |85[ ) and 
(|86[) are in qualitative agreement since ~ g 2 T 2 and both rates are of order e 2 g 2 T 4 /q . In fact the 
enhancement mechanism discussed above is realized in M, |8[ |9| by the factor T 2 /M^. We interpret this 
pole in as a use of the approximation of Eq. (]83| ) , extended to the region where L ~ T in order to 
get the logarithmic factor ln(T 2 /M g ). This leads to a very different limiting behavior when — ► 
since in our approach we recover a logarithmic singularity while in the semi-classical approach one 



obtains a pole in |J. Another difference between Eq. (85) and the results of Cleymans et al. is 
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that the transverse gluon exchange gives a finite contribution which is equally important as that of the 
screened longitudinal one considered in |7|, ||, g]. As a consequence, the use of the static approximation 
to describe the interaction of the quark in the plasma is not justified in thermal field theory. 



6 The slightly virtual case: g 1 < Q 2 /ql <C 1 



We now consider the case (ii) of Sec. 3.5 where the virtuality Q 2 is non negligeable but remains small. 
This can be treated in a simple way thanks to the introduction of an effective fermion mass. 



6.1 Expression of ImlP ^ 

As already mentioned, in this case we still have a clear hierarchy in the various powers of Q 2 lq\ 



< 1. 



The consequence of this hierarchy is that the dominant term of the numerator remains the same as 
in Eqs. (]45|) , and that it can still be approximated by Eq. (^). In what concerns the denominators, 
we saw in Eqs. (|38|), (|39| ) that their expression remains unchanged provided one replaces the fermion 
thermal mass by an effective mass taking into account the effect of the virtuality of the emitted 
photon M 2 S = + Q 2 r 2 /qQ. Even before writing the expression of the imaginary part of the photon 
self-energy, we can notice that the positive virtuality Q 2 > has the consequence of increasing the 
effect of the fermion thermal mass. Another point to be looked at is the influence of the constraint 
cos 9' £ [—1,1] on the phase space. By examining the denominators entering in the rate, we easily 
see that the relevant values of u in the integration over u are of order of u* ~ M 2 s /r 2 . As far as 
Q 2 /qq <C 1, this collinear regulator is much smaller than 1. Therefore, we can still approximate 
cos#' ~ lo/l. Therefore, the effect of the constraint cos6' £ [—1,1] is still to confine L 2 to negative 
values. Finally, in complete analogy with Eq. the imaginary part of the self-energy of a virtual 
photon reads: 



Imn%(Q) «(-l) 



eV 1 
T 87T* q 



+1 



2\2 



dm' (r) I dl I dx n B {lx) p T L (l,lx)(l — x ) 



u + 



M 2 , 
cf 

W 2 



du 



u + ^pr^ — j F 



the notations being the same as in Eq. (47) 



1/2 ' 



(87) 



6.2 Reduction of the expression 



The simplification of this imaginary part can be performed by the same tools as in Sec. 4.2. In 
particular, since the relevant values of the angular variable u are much smaller than one, the same 
approximations are valid. Despite that similarity, an essential difference lies in the fact that the 
integration over r cannot be factorized since the effective fermion mass now depends on r. If one 
introduces the dimensionless variables: 

V = f 
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— L 2 

w 



M? s (v) 



ImU 

I T r{X,V) = 



~ ReII Ti (x) 



Eq. ( p7|) takes the form 



+00 

1 



r*/T 



e« + 1 



1 to*(t;) 

I x T ' LK J (w + R TL (x,v)) 2 + (I TL (x,v)) 2 







where, compared to Eq. (|55|), the integration over v cannot be factorized because of the v dependence 
contained in the quantities w*, R T L and I TL . 

6.3 Discussion of the result 

When the virtuality Q 2 > is switched on, the quantity which is worth studying is now the integral: 



+00 1 



6 f 1 

T = —K / dVV 2 



1 



e v + 1 



r*/T 
w* (v) 



x / dw vWC^tarm-lyyR+4) ( } 

where the explicit prefactor 6/-7T 2 is introduced to compare this definition for J TL with the previous 
one in the limit where Q 2 = and where the integral over v can be factorized. 

A priori, the functions J T L depend of the four dimensionless quantities r*/T, Z*/M oc , Moo/nig and 
Q 2 T 2 /qQm 2 . Nevertheless, since the weight function p(v) = v 2 e v /(e v + l) 2 is peaked around v ~ 1, 
the small values of v will be negligeable in the result, which means that in fact this result does not 
depend on the value of the cut-off r*, which thus can be taken to zero. 

Another question which can be answered very simply by using the study done for the production 
of real photons concerns the absence of IR divergences when the emitted photon is virtual. Indeed, 
we saw earlier that the IR finiteness was in fact due to the presence of a fermion thermal mass. In 
the present case, the effect of this fermion mass is enhanced by the presence of a positive virtuality 
Q 2 > 0, since M cff > M^. Therefore, it is obvious that the production rate of virtual photons will be 
IR safe, even for the contribution of the transverse gluon, and without the need of a magnetic mass. 

To end with generalities, one can notice that the functions J T L are decreasing functions of the 
virtuality Q 2 . Indeed, if Q 2 increases, the effective mass M c g(r) increases uniformly in r. Then, 
according to the study performed in the case of real photons, the functions J T L should decrease. 
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Concerning the dependence on the UV cut-off I*, one can evaluate the contribution of the region 
I 6 [I*, +oo[ as being: 



M cS (v) 



hOO 



Neglecting the logarithms, the contribution of the hard region [I*, +oo[ is of order (m g /l*) 2 , and must 
be compared to J TL in order to know if we can make I* = +00. If we recall that the values which 
contribute dominantly in the integral over v are of order 1, we see that the influence of the virtuality 
Q 2 > on J T L is controlled by the comparison of and Q 2 T 2 /qfi (see Eq. (|40|)). The results are 
as follows: 

(i) If Q 2 /ql) <C M^/T 2 ~ g 2 , the virtuality Q 2 has no effect and can be neglected: this case 
corresponds to the previous section. The functions J T L are of order 1, and the contribution of 
[I*, +00 [ is negligeable as far as I* 3> m g . It means that the hard region around I ~ T does not 
contribute. 

(ii) If Q 2 lq\ ~ g 2 , we are in the intermediate region where the virtuality of the radiated photon and 
the fermion thermal mass have comparable effects. The functions J TL are still of order 1, and 
the conclusion of the previous case remains unchanged. 

(hi) If g 2 -C Q 2 /qQ, the effect of the virtuality becomes much more important than the effect of 
the fermion thermal mass M^, so that M 2 S ~ Q 2 r 2 /q 2 ,. Moreover, in that case, we have 
nig <C M e g(v ~ 1). By using the same tools as in the previous section, one obtains easily the 
simplified expression: 
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3 f . 2 1 
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e v + 1 




Xx,v) (92) 



and, separating the transverse and longitudinal terms, one gets: 



TV ; 87r Q 2 T 2 \qlmlj 

J L (l* -»■ 00) » -2J T (l* 00), (93) 



where the i? dependence is neglected inside the logarithm. We have checked that this approxi- 
mation is in agreement with numerical estimates of Eq. (^) . We are now in position to compare 
the contribution from [Z*,+oo[ with that of [0,/*]. To make the discussion more precise, let us 
assume that Q 2 /q^ = g a and I* = g^T, with a, j3 > 0. Then, neglecting the logarithmic factors 
in our comparison, we have: 

J T,L « J T,L 
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J' TL « J T Jl* ^oo) 

g a « g w 

(3 < a/2 since j< 1. 



(94) 



This means that for fixed Q 2 /q 2 = j a < 1, the dependence on /* disappears as long as I* > g@T 
with < P < a/2, and so the contribution to the result of the region beyond such a I* is 
negligeable. Therefore, in the region which contributes dominantly to the result, we still have 
since j3 > 0. This means that our kinematical approximations are safe. Of 
course, this works since we had Q 2 jq 2 . <C 1 so that a is strictly positive. When the ratio Q 2 jq 2 . 
approaches the value 1, the negligeable region around I ~ T becomes smaller and smaller, so 
that the precision of our approximations decreases. Another way to say the same thing is that 
given a fixed value of g <C 1 and given an accuracy one wants to ensure, there exists an a > so 
that Q 2 jq 2 . must remain lower than g a . In order to reach values of Q 2 jq 2 . closer to 1, one has to 
relax the constraint on the accuracy or to consider a smaller g. We emphasize that, although we 
have introduced here fractional powers of the coupling constant g in order to quantify when our 
approximations for evaluating the bremsstrahlung diagrams are valid, this does not mean that 
the perturbation expansion, in integer powers of g, breaks-down; nevertheless, before calculating 
next-to-leading order diagrams, one should improve the approximations for the leading ones. 

A general consequence of this analysis is the fact that the relevant scales in our problem are 
not only gT and T like in the standard HTL framework, but also the scale of the effective mass 
M cS (v ~ 1) ~ + Q 2 T 2 jq 2 which is intermediate between gT and T. This is a consequence of the 
fact that the regularization of all our divergences, including the infrared one, is done by this effective 
mass rather than by the resummation of the HTLs in the gluon propagator. 

The case Q 2 jq\ ~ 1 does not fit in these approximations since in that case a is of order zero, so 
that (3 is also of order zero and even the hard values of I contributes to the result. 

To summarize this section, the functions J T L have been plotted as a function of log(Q 2 /<7o) on 
Fig. ^ over the whole range of admissible values of Q 2 . We observe a clear flattening of the curves 
for small enough Q 2 /q 2 , indicating that the dependence on Q 2 becomes negligeable. In the opposite 
direction, we see that the enhancement of the rate decreases if Q 2 — > q 2 . The enhancement (ratio 
between the values of J TL at Q 2 /q 2 = and at Q 2 jq\ ~ 1) decreases also if the coupling constant 
g increases. We note also that the longitudinal and transverse contributions are almost equal in 
magnitude. 



6.4 Reduction using sum rules 



Again, it is possible to reduce Eq. ( p7| ) by performing exactly the integral over x thanks to sum rules. 
Since now we cannot factorize the integral over r, the result given by sum rules is a two-dimensional 
integral. The techniques to achieve that are exactly the same as in Sec. 4.3, and therefore will not be 
reproduced here. We still need to distinguish three cases. For example, for the phenomenologically 
relevant case, 3/8 < M^/m 2 the quantities 1 - 4M c 2 ff (r)/ReII Ti (x), 1 - 4M 2 s (r)/(uj 2 - I 2 ) and 
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1 — 4M 2 ff (r)/m 2 are negative for every value of r. The integration over u' gives in the transverse case: 
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For the other cases 1/4 < M^/m 2 < 3/8 and M^/m 2 < 1/4 similar calculations can be carried out 
but they are much more complex since the signs which control the expressions depend on both v and 
I. 



7 The case: Q 2 /ql - 1 

This case will not be studied in detail in this paper; we discuss only why it is very different from the 
previous ones, and what is the expected order of magnitude of the production rate. This case differs 
from the previous two in the following respects: 

(i) Since Q 2 /qfi ~ 1, all powers of this ratio should now be kept in the expansion of the quantities 
we need. 

(ii) The relevant values of the angular variable u are now of order u* ~ 1. Therefore, the collinear 
approximation (which means keeping only the term in vP in the numerator) cannot be used. 
For the same reason, the approximation cos#' ~ Iq/1 is no longer valid. This means that the 
contribution of the domain I? > is as important as the contribution of the Landau damping 
L 2 < 0. Physically, this corresponds to processes involving the Compton effect (a thermalized 
on-shell gluon is absorbed by the moving quark, which emits a photon). 

(iii) The scattering of the quark is now sensitive to parton exchanges of hard momentum. As a 
consequence, the approximation Iq,1 <C T cannot be used in this case. Moreover, since our 
argument to neglect the diagram of Fig. Q-(c) was essentially based on the possibility to neglect 
a Fermi-Dirac statistical weight in front of a Bose-Einstein one when Iq T, it is now impossible 
to apply it. Therefore, this diagram contributes now at the same order. 

Extrapolations of the previous case (despite its incompleteness) based on Eqs. (|93"|), as well as 
preliminary estimations of the order of magnitude for the production rate of a static photon (Q 2 /q® = 
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1), seems to indicate that the bremsstrahlung gives a contribution Imll 7 ^ ~ e 2 g 4 T 3 /qo- Therefore, the 
diagram of Fig. [j]-(a), already considered by Braaten, Pisarski and Yuan, is not the only contribution 
to the production of virtual photons. 



8 Summary 

We have studied the production of photons of energy qo <C T by a plasma. For virtualities verifying 
Q 2 /Qq ^ 1) this photon emission is dominated by the bremsstrahlung process. In this calculation 
of the photon production rate, the HTL framework appeared to be insufficient in two respects: first 
of all, the standard HTL power counting breaks down in extracting the dominant contributions to 
this rate. Indeed, the diagram we considered is supposed to be sub-dominant in the HTL expansion, 
but turns out to dominate over the diagrams considered in the HTL framework. The reason for 
this failure lies in a strong enhancement of our diagrams due to collinear sensitivity. Secondly, once 
the relevant contributions have been isolated, the HTL resummation of soft lines is insufficient to 
regularize collinear divergences. This problem is solved by the application of a recent extension of 
the HTL resummation program. The imaginary part of the polarization tensor we get, proportional 
to the production rate, goes from Imn^ ~ e 2 g 2 T 3 /qo when Q 2 /qQ <C g 2 to Im ~ e 2 g 4 T 3 /qo 
when Q 2 /qQ ~ 1. At the same time, when the ratio Q 2 jq 2 . increases, the bremsstrahlung process 
becomes more and more sensitive to gluon exchanges of momentum intermediate between the soft 
scale gT and the hard scale T, which is again an important qualitative difference with the standard 
HTL calculations. 

When the virtuality Q 2 /q 2 , approaches 1, the enhancement disappears, but the bremsstrahlung 
process remains of an order comparable to the contributions already calculated by Braaten, Pisarski 
and Yuan. This may be interpreted as follows: the rate obtained within the HTL expansion by BPY 
is unexpectedly subleading due to a partial cancelation when taking the trace of the polarization 
tensor H^ u . Therefore, other terms subleading in the HTL power counting rules may contribute as 
well. 
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10 Appendix 

To discuss the asymptotic behavior of the integrals J TL defined in Eq. (|70|), it is convenient to split 
the integral over w in J T L in order to obtain: 
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In each of the two terms, we will use the small w approximation: 



(96) 



w/(w + 4)tanh 1 ^Jw/(w + 4) p^ ^ 



1/(1 + 4«;)tanh-y 1/(1 + 4«;) « (97) 

Moreover, since the dominant behavior of the integrals J T L is controlled by the small x region, we 
will use the following approximations 

R T {x) « — £ I.faO ~ V 

~ , , 3m1x 2 ~ , 3m1x , „ 

R A X ) ~ *,o It(?) ra TT^T- (98) 

tV ; x«i 2M^ tV ' x<l 4M^ v ; 

The only supplementary ingredient we need is to recall that the two functions R T L (x) and I T L (x) are 
proportional to the ratio m^/M^, which determines their order of magnitude, compared to w which 
is of order 1. After a bit of algebra, one obtains the Eqs. (f7l|). 
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Figure Captions 

Fig. |l] Contributions to the soft photon production rate with soft internal lines. 

Fig. |2] a,b,c: matrix elements contributing to the imaginary part of diagrams (a) and (c) of the previous 
figure; d,e: matrix elements contributing to the imaginary part of the gluon tadpole diagram in 
Fig. § 

Fig. |H Contributions to the soft photon production rate with hard internal fermion lines - (a): vertex 
insertion; (b): self-energy insertion. 

Fig. | l-4M^/Ren Tji (x). 

Fig. |5] Effect of the magnetic mass on J T . 

Fig. |H J TL as a function of log((5 2 /<7o)- 
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Figure 2: (a),(b),(c): matrix elements contributing to the imaginary part of diagrams (a) and (c) in 
Fig. |]; (d),(e): matrix elements contributing to the imaginary part of the gluon tadpole diagram in 
Fig. |l[ The solid dots indicate effective propagators or vertices. The emitted photon has momentum 
much less than T. The other lines carry hard momentum. 
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(b) 




/' 



Figure 3: Contributions to the soft photon production rate with hard internal fermion lines - (a): 
vertex insertion; (b): self-energy insertion. The gluon propagator is an effective one. The fermion in 
the loop is hard. 
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Figure 4: Solid line: l-4M^/ReII T (:r). Dashed line: l-4M^,/ReU L (x). Dotted line: l-AM^/mj. 
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Figure 5: Effect of the magnetic mass on J T . Solid line: (M^/m^) 2 = 1. Dotted line: (M 00 /m l 
0.1. 
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Figure 6: J T L as a function of log(Q 2 /qfi), with {m^/M^) 2 fixed to 1.5 (i.e. N = 3 colors and N f = 3 
light flavors). Solid curves: (m g /T) 2 = 0.005 or g = 0.1. Dashed curves: (m g /T) 2 = 0.1 or 5 = 0.44. 
Dotted curves: (m g /T) 2 = 2 or 5 = 2. 
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